Effective conductivity of composites of graded spherical particles 
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Abstract 

We have employed the first-principles approach to compute the effective re- 
sponse of composites of graded spherical particles of arbitrary conductivity 
profiles. We solve the boundary- value problem for the polarizability of the 
graded particles and obtain the dipole moment as well as the multipole mo- 
ments. We provide a rigorous proof of an ad hoc approximate method based 
on the differential effective multipole moment approximation (DEMMA) in 
which the differential effective dipole approximation (DEDA) is a special case. 
The method will be applied to an exactly solvable graded profile. We show 
that DEDA and DEMMA are indeed exact for graded spherical particles. 

PACS Number(s): 77.22.-d, 77.84.Lf, 42.79.Ry, 41.20.Cv 



Typeset using REVTeX 



* Corresponding author. Electronic mail: kwyu@phy.cuhk.edu. hk 



1 



I. INTRODUCTION 



In functionally graded materials (FGM), the materials properties can vary spatially 
These materials have received considerable attention in various engineering applications [1] . 
The variation in the composition yields material and microstructure gradients, and makes the 
FGM very different in behavior from the homogeneous materials and conventional composite 
materials [1,2]. The great advantage is that one can tailor the materials properties via the 
design of the gradients. Over the past few years, there have been a number of attempts, both 
theoretical and experimental, to study the responses of FGM to mechanical, thermal, and 
electric loads and for different microstructure in various systems [1-8]. In Nature, graded 
morphogen profiles can exist in a cell layer [4]. In experiments, graded structure may 
be produced by using various approaches, such as a three-dimensional X-ray microscopy 
technique [7], deformation under large sliding loads [8], and adsorbate-substrate atomic 
exchange during growth [6]. Thus, gradation profiles exist in both natural materials and 
artificial FGM. Interestingly, gradation profiles may further be controlled according to our 
purpose, such as a specific power-law gradation profile and so on. It has reported recently 
that the control of a compatibility factor can facilitate the engineering of FGM [5] . 

There have been various different attempts to treat the composite materials of homo- 
geneous inclusions [9] as well as multi-shell inclusions [10-13]. These established theories 
for homogeneous inclusions, however cannot be applied to graded inclusions. To this end, 
we have recently developed a first-principles approach for calculating the effective response 
of dilute composites of graded cylindrical inclusions [14] as well as graded spherical parti- 
cles [15]. The electrostatic boundary- value problem of a graded spherical particle has been 
solved for some specific graded profiles to obtain exact analytic results. Along this line, 
exact analytic results are available so far for the power-law graded profile [14,15], linear 
profile [14], exponential profile [16,17] as well as some combination of the above profiles [18]. 
For arbitrary graded profiles, we have developed a differential effective dipole approxima- 
tion (DEDA) to estimate the effective response of graded composites of spherical particles 
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numerically. The DEDA results were shown to be in excellent agreement with the exact 
analytic results [15]. However, the excellent agreement is difficult to understand because 
DEDA method was based on an ad hoc approximation [19]. 

The object of the present investigation is two-fold. Firstly, we will extend the first- 
principles approach slightly to deal with graded particles of arbitrary profiles. We will solve 
the boundary-value problem for the polarizability of the graded particles and obtain the 
dipole moment as well as the multipole moments. Secondly, we will provide a rigorous proof 
of the ad hoc approximation from first-principles. To this end, we extend the proof to mul- 
tipole polarizability and derive the differential effective multipole moment approximation. 
As an illustration of the method, application to the power-law profile will be made. Thus, 
both DEDA and DEMMA are indeed exact for graded spherical particles. 

II. FIRST-PRINCIPLES APPROACH 

In this work, we will focus on a model of a graded conducting particle, in which the 
conductivity of the particle varies continuously along the radius of the spherical particle. We 
consider the electrostatic boundary-value problem of a graded spherical medium of radius a 

subjected to a uniform electric field E applied along the z-axis. For conductivity properties, 

— * — » — * — * 
the constitutive relations read J = (Ti(r)E and J = a m E respectively in the graded spherical 

medium and the host medium, where o~i(r) is the conductivity profile of the graded spherical 

medium and a m is the conductivity of the host medium. The Maxwell's equations read 

V-J = 0, Vx£ = 0. 

To this end, E can be written as the gradient of a scalar potential $, E — — V$, leading to 
a partial differential equation: 

V • [cr(r)V$] = 0, (1) 

where a(r) is the dimensionless dielectric profile, while a(r) = (Ji{r)/a m in the inclusion, 
and o~(r) = 1 in the host medium. 
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In spherical coordinates, the electric potential $ satisfies 

1 d ( 2 . 1 d ( . n . 1 d ( . n 

^ ( r ^ aFj + r^Tem l^Ww) + r^Ted^ \ a{r) d^) = °" (2) 

As the external field is applied along the z-axis, $ is independent of the azimuthal angle 
(p. If we write $ = /(r)0(0) to achieve separation of variables, we obtain two distinct 
ordinary differential equations. For the radial function /(r), 

Jr(^(0^)-i(i + iMO/(0 = o, (3) 

where / is an integer, while Q(9) satisfies the Legendre equation [9]. Eq.(3) is a homogeneous 
second-order differential equation; it admits two possible solutions: ff(r) and ff(r) being 
regular at the origin and infinity respectively. Exact analytic results can be obtained for a 
power-law profile [14,15], linear profile [14], and exponential profile [16]. The general solution 
for the potential in the spherical medium is thus given by 

oo 

Mr,0) = £[AJ+(r) +B i / i -(r)]P J (cos0). (4) 

1=0 

In the host medium, the potential is given by 

oo 

®m(r,9) =£[<V + Ar- (I+1) ]Pj(cos0). (5) 

1=0 

Thus the problem can be solved by matching the boundary conditions at the spherical 
surface. 

III. BOUNDARY- VALUE PROBLEM 

By virtue of the regularity of the solution at r = 0, the general solution inside the particle 
becomes: 

oo 

$,(r,0) = X>/* + W(cos0), r<a, (6) 



1=0 



The external potential (r > a) is: 
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$ m (r,9) = -E r cos 9 + J2 D i r ~ (l+1)p i( cos0 ), r > «• (7) 

1=0 

We can rewrite rcos9 as J2i Snr 1 Pi(cos 9) . Thus Ci = —E Sn and no multipole moment 
will be induced in a uniform applied field. However, in a nonuniform applied field as in di- 
electrophoresis of graded particles, multipole moments will be induced. Matching boundary 
conditions at r = a, 

$i(a) = $ m (a), <r(a)$;(a) = $^(a), 

where prime denotes derivatives with respect to r, we obtain a set of simultaneous linear 
equations for the coefficients 

A l f+(a) = -E a l 5 ll + D l a-« +1 \ (8) 
<7(a)A,/+'(a) = -Eola^Sn - D t (l + l)a^ l+2 \ (9) 

Solving these equations, 

(2/ + 1 Vflofa (21 + 1 Vflpfa 

' aa(a)/+'(a) + (I + l)tf-(a) /+(a)[/(F z + 1) + 1] ' 1 j 

a 2;+1 [aa(a)/ ; +/ (a) - l/, + (a)]£ fti _ a 2 ' +1 Z(F, - ljifrfa ( . 

aa{a)f+'{a) + {l + l)f+{a) l(F t + 1) + 1 ' 1 } 

where 

cri(o) a/; +/ (a) 
°m lji{a) 

When the radial equation is solved for a specified graded profile, the potential distribution 
can generally be expressed in terms of ff(r) and F\. For a uniform applied field, the potential 
becomes 



3aA + (r)ff 
/+(a)(F 1 + 2) 



*i(r,0) = - o^ cos ^ r <«> ( 13 ) 



$ m (r, 0) = -F r cos + cos 0, r > a. (14) 

The second term in the potential in Eq.(14) can be interpreted as the potential due to an 
induced dipole moment p = cr m bia 3 E . Thus we identify the dipole factor: 



5 



h F1 ~ 1 J? ^( Q ) a fl'( a ) n ^ 

For a homogeneous (non-graded) particle, the well known result recovers: 

h = 



Oi + 2(7 m 

Thus, F 1 can be interpreted as the equivalent conductivity ratio of the graded spherical 
particle. 



IV. DIFFERENTIAL EFFECTIVE MULTIPOLE MOMENT APPROXIMATION 

We should remark that in general Fi is the equivalent conductivity ratio of the /th 

multipole moment. Let us extend the definition [Eq.(12)] slightly for a graded spherical 
particle of variable radius r: 

W = *Mgff (16) 

Thus, the multipole factor reads: 

W(r) - 1) 

= /(fiM + D + i- (17) 

Physically it means that we construct the graded particle by a multi-shell procedure [19]: we 
start out with a graded particle of radius r and keep on adding conducting shell gradually. 
The change in F\ and bi can be assessed. In this regard, it is instructive to derive differential 
equations for F;(r) and 6/(r). Let us consider (ignoring the superscript + and subscript / in 
f^ij-) for simplicity): 



2 , ./'I'') 



fir) 



1 d r 2 , w„ m r 2 a(r)f'(r) 2 
[r 2 a(r)f(r)} 



lf(r) dr lf{r)' 



From Eq.(3) and Eq.(16), we obtain the differential equation: 

d IF(r) 2 

T [rF l (r)] = (l + l)v(r)--!$-L, (18) 
ar a(r) 

which is just the generalized Tartar formula [20]. From Eq.(17) and Eq.(18), we obtain the 

DEMMA: 
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Mr) 



1 



1 + b l (r)l) + (b l (r)-l)la(r)] 



dr (21 + l)m(r) 

x [(6,(r) + / + 6,(r)Z) - (6,(r) - 1)(Z + l)<r(r)]. 

When Z = 1, we recover the Tartar formula and DEDA [19]: 

-[^(01=2.(0-^-, 



(19) 



(20) 



^i(r) 
dr 



■[(26 1 (r) + l) + (6 1 (r)-l)<r(r)] 



3rcr(r) 

x[(26i(r) + l)-2(&i(r)-l)<7(r)]. 



(21) 



Let us consider a graded particle in which the conductivity profile has a power-law 

dependence on the radius, a(r) = cr k , with k > where < r < a. Then the radial 
equation becomes 

d 2 / fc + 2 d/ Z(Z + 1) 



dr 2 ^ r dr 



-f = o. 



(22) 



As Eq.(22) is a homogeneous equation, it admits a power- law solution [15], 



f(r) = r s . 



(23) 



Substituting it into Eq.(22), we obtain the equation s 2 + s(k + 1) — 1(1 + 1) = and the 
solution is 



4(Z) = - -(fc+l)±^/(fc + l)2 + 4Z(Z + l) 



(24) 



There are two possible solutions: 



ft(r)=r*W, /f(r)=r 



**(0 



Thus the Z-th order equivalent conductivity becomes: 

^(a) q/Qo) s k + (l) k 

When k — > 0, s+(Z) — ► Z, F; — > c, the result for a homogeneous sphere recovers. 



(25) 



V. DISCUSSION AND CONCLUSION 



Here a few comments are in order. We have employed the first-principles approach to 
compute the multipole polarizability of graded spherical particles of arbitrary conductivity 
profiles and provided a rigorous proof of the differential effective dipole approximation. 

We are now in a position to propose some applications of the present method. We may 
attempt the similar calculation of the multipole response of a graded matallic sphere in 
the nonuniform field of an oscillating point dipole at optical frequency. The graded Drude 
dielectric function will be adopted [21]. The similar approach may also be extended to 
anisotropic medium with different radial and tangential conductivities [22]. Similar work 
can be extended to ac electrokinetics of graded cells [23] . We can also study the interparticle 
force between graded particles [24]. 

In summary, we have solved the boundary-value problem for the polarizability of the 
graded particles and obtained the dipole moment as well as the multipole moments. We 
provided a rigorous proof of the differential effective multipole moment approximation. We 
showed that DEDA and DEMMA are indeed exact for graded spherical particles. Note that 
an exact solution is very few in composite research and to have one yields much insight. Such 
solutions should be useful as benchmarks. Finally, we should remark that the exact deriva- 
tion of DEDA and DEMMA is for graded spherical particles only. For graded nonspherical 
particles, these ad hoc approaches may only be approximate. 
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